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Graphene on metals 

1. Introduction 

Graphene (gr), a two-dimensional layer of carbon 
atoms arranged in a honeycomb lattice [Fig. [lja)], 
demonstrates a variety of unique electronic and 
transport properties. The discovery of such fascinating 
phenomena as very high electron and hole mobilities, 
ambipolar electric field effect, integer and half-interger 
quantum Hall effects for electrons and holes, etc. [l][2] 
attracts a lot of attention from the different fields of 
solid state physics and chemistry. 

The hexagonal lattice of graphene consists of two 
equivalent carbon sublattices [A and B in Kig. [lja)]. 
Carbon atoms in this structure are sp 2 hybridzed (one 
2s and two 2 p orbitals), that leads to the formation 
of the planar structure with strong a bonds and the 
distance of 1.42 A between carbon atoms. Completely 
filled a bands for the infinite graphene layer forms deep 
valence band levels of graphene. The 2 p z orbitals of the 
neighbouring carbon atoms are perpendicular to the 
plane of graphene and their overlap above and below 
this plane leads to the formation of the 7r valence band. 
This band is half filled and the unique property of the 
electronic structure of graphene is that the n and 7r* 
bands touch at the corners of the hexagonal Brillouine 
zone of graphene (at the K points) directly at the 
Fermi level (Ep) [Fig. [ljb,c)]. The band dispersion, 
E(k ), of the 7r states in the close vicinity of Ep is linear 
and described via the Dirac equation E = Hvpk, where 
vf ~ 1 x 10 6 m/s is the Fermi velocity (“speed of light” 
for the massless Dirac fermions). The point where 7r 
bands intersect at I\ is called a Dirac point. Thus, 
graphene is a semimetal in the undoped case. Due 
to the low density of the valence band states (DOS) 
of graphene around Ep and the linear dependence of 
DOS on energy, the electronic structure of graphene 
can be strongly modified in a controllable way via 
attaching different species to graphene that can lead 
to different technological applications of graphene, like 
touch screens or gas sensors [3}j5]. 

First transport experiments on graphene were 
performed on the single graphene flakes exfoliated 
through the so-called “scotch-tape” method from 
bulk graphite. However, despite the high crystalline 
order and transport mobility of these flakes, they 
cannot be used in the technological process due 
to their small size and uncontrollable preparation 
procedure. In a series of experiments it was 
shown that graphene layers of very high quality can 
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Figure 1. (a) Crystallographic structure of graphene with 

two equivalent carbon atoms, A and B, in the unit cell 
marked by the solid-line black rhombus, (b) Two-dimensional 
electronic struc ture of graphene [E(k x , k y )\ in the tight-binding 
approximation [l 0 . (c) The DFT calculated band structure of 
graphene in the vicinity of Ep along the main directions in the 
hexagonal Brilloine zone (shown as an inset) El- The dispersion 
of 7r and 7r* states is highlighted by the thick line. 

be prepared on SiC semiconducting substrates ||. 
These graphene layers on SiC demonstrate transport 
properties comparable to those for the exfoliated 
graphene |7|[8] and it was shown that the electronic 
properties of graphene/SiC interface can be tailored 
in a wide range of properties via adsorption or 
intercalation of different species |6 ; . Later it was 
demonstrated that graphene nanoribbons can be 
prepared on the templated SiC substrate using scalable 
photolithography and microelectronics processing and 
this technology leads to the fabrication of the graphene 
transistors with density of 40,000/cm 2 |9 . 

However, as was shown in the recent works the 
most promising method for the preparation of the 
graphene layers, which can be further used in the 
technological applications, is the synthesis of single- 
and multilayer graphene on the metallic support 
via chemical vapour deposition (CVD). This method 
combined with roll-to-roll method allows production 
of the graphene layers on metallic substrate and then 
transfer on the polymer support up to meter sizes [4]. 
In the end of 2009 the full technological process was 
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presented that leads to the fabrication of the first 
graphene-based touch screens. Later, the optimisation 
of this technology allowed to fabricate the first mobile 
phone with the screen fully made on the basis of 
graphene These graphene layers demonstrate the 
electron mobility of 5100 — 5200 cm 2 V _1 s _1 at room 
temperature, visible light transparency of 97.75% at 
550 ran, and quantum Hall effect [4}[5]. 

Graphene (known as monolayer of graphite in a 
former time) overlayers on metallic surface are in focus 
of surface science research for many years starting 
from middle of the 60s [l2 


graphene layer on Pt(100) 


18 


12 


, when the multidomain 
and the moire graphene 


structure on Ru(0001) [15] were identified in the low- 
energy electron diffraction (LEED) experiments. At 
that time the main interest to this topic was motivated 
by the studies of the catalytic activity of the clean 
d-metal surfaces and in most cases graphene was 
considered as a poison layer blocking their reactivity. 

The discovery of the fascinating transport proper¬ 
ties of graphene renewed interest to the investigation 
of the graphene/metal interfaces. Here the studies are 
ranged from the fundamental problems, like the correct 
description of the relatively weak interaction between 
graphene and metal surface, which in most cases leads 
to the drastic modification of the electronic structure of 
graphene, to the more practical issues, like the prepa¬ 
ration of the ordered arrays of clusters on the graphene 
moire on the surface of 4 d and 5 d metals, which then 
can be used as a model system for the studies of e. g. 
catalytic properties of these systems. 

The present manuscript gives a timing review 
of the recent experimental and theoretical results on 
the graphene/metal interfaces. Starting from the 
consideration of the methods of graphene’s growth 
on metallic surfaces, it follows with the overview of 
the crystallographic and electronic properties of these 
systems, including the experiments on the modification 
of graphene/metal via adsorption and intercalation 
of different species and formation of graphene nano¬ 
objects. This review includes several examples of 
the representative results obtained with the main 
surface science methods (STM, AFM, LEEM-PEEM, 
NEXAFS, XPS, ARPES) and tries to make links 
between them with the help of the state-of-the-art 
density functional theory (DFT) calculations. 


2. Experimental methods 

Graphene layers on metallic surfaces are pure 2D 
objects and they are ideal objects for the application 
of different surface science techniques. Low-energy 
electron microscopy (LEEM) and scanning probe 
microscopy (SPM) methods give information about 
morphology of the system and the electronic structure 


on the local scale (from pm to atomic scale). 
The photoelectron spectroscopy based techniques 
(NEXAFS and PES) allow to get information about 
the electronic structure on the macroscopic scale (from 
several hundreds nrn to mm scale). This section gives 
a short overview of the listed techniques pointing their 
(dis-) ad vant ages. 

2.1. LEEM/PEEM 

LEEM is a parallel imaging technique, which uses 
elastically backscattered electrons to image crystalline 
surfaces and interfaces [Fig. [2]a)]. Due to the large 
electron backscattering cross section of most materials 
the intensities of the diffracted and/or reflected beams 
are quite high, that makes LEEM the ideal technique to 
image in video rate dynamic processes such as surface 
reconstructions, epitaxial growth, step dynamics, self¬ 
organization and others. LEEM enables to reach a 
lateral resolution of less than 10 nm. Most importantly, 
it offers several structure sensitive complementary 
imaging and diffraction methods to probe the surface. 
The detailed description of the operation principles of 
this method can be found elsewhere [l9] - [22] . 

A beam of high energy electrons (10 — 20keV), 
decelerated by the retarding potential of the objective 
lens, impinges in normal incidence on the surface, 
with energy in the range 0 to few hundreds eV. The 
beam energy is varied by adjusting a bias voltage 
between sample and the emitter. The elastically 
backscattered electrons are then reaccelerated through 
the objective lens, following the inverse pathway. The 
objective produces a magnified image of the surface 
in the beam separator, which is further magnified by 
several additional lenses in the imaging column of the 
instrument. This image is projected onto an imaging 
detector with microchannel plate and phosphorous 
screen, and finally acquired by a CCD camera. Along 
with the real space imaging, a LEEM microscope is 
also capable of the reciprocal space imaging. 

Depending on the diffracted beam selection either 
bright field of dark field imaging in LEEM is 
possible. If the zero-order diffracted beam (“00” 
beam) is selected by a contrast aperture positioned 
in the diffraction plane then the bright-held LEEM is 
performed (defocusing of the beam allows to convert 
the phase difference, appeared due to the reflection 
from different heights, into an amplitude difference 
that gives a possibility to image steps at the surface). 
In this case the pure structural contrast is obtained. In 
case of overlayers on the surface, the interference of the 
backscattered beams from the surface and interface can 
produce the maxima and minima in the intensity that 
allows to measure a local thickness of the overlayer. If 
the secondary diffracted beam is selected then a dark- 
held image of the surface is produced. Here all areas 
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Figure 2. Principal schemes of (a) low-energy electron microscopy (figure is taken from Ref. 19] with permission), (b) scanning 
probe microscopy/spectroscopy (STM/AFM), and (c) photoelectron spectroscopy (XPS/ARPESj and near-edge x-ray absorption 
spectroscopy (NEXAFS/XMCD) experiments. 


that contribute to the formation of the selected beam 
appear bright. 

In mirror electron microscopy (MEM) the surface 
is illuminated with electrons at very low energy, so 
that the electrons interact very weakly with the surface 
(this occurs at the transition MEM-LEEM). Under 
these conditions the contrast is due to work function 
differences and topography variations. MEM allows 
non-crystalline samples to be imaged. 

In case if photon source is used instead of the 
electron gun, a LEEM can be used for the imaging of 
photoemitted electrons. If the entrance slit is placed on 
the entrance of the the prism array (MPA) on the way 
of the emitted electron than the use of the reciprocal 
imaging produce the fc-space resolved photoelectron 
intensity map allowing to perform local imaging of the 
electronic structure. 


2.2. Scanning probe microscopy methods 


SPM methods combined with the corresponding 
spectroscopy add-ons give an information about 
morphology of the studied objects and the electronic 
structure down to the atomic scale [Fig. §b)] [2§ 


25 


In scanning tunnelling microscopy (STM) the 
feedback is regulated via setting the value of the 
tip position above the surface, z, that regulates the 
tunnelling current, It, which exponentially depends 
on that gives the distribution of the local density of 
states at the surface, p(x,y,z). Here Ip(x,y, z,V) oc 
f p(x,y,z,E)dE , where the integration is performed 
in the energy range between the Fermi levels (Ep) of 
tip and sample, when they differ by the value eUp, 
where Up is the bias voltage. One can perform the 
differentiation of the I(x,y, z,Up) signal with respect 
to the bias voltage (scanning tunnelling spectroscopy, 


STS). This can be carried out with lock-in technique 
and gives a direct information about the local density 
of states. If such measurements performed in the 
scanning mode at different bias voltages, then it allows 
to observe the so-called electron density standing waves 
of different periodicities that can be used to obtain 
the electron dispersion relation, E(k), of the surface 


electronic states 26 -28 


In the modern atomic force microscopy (AFM) 
experiments, the oscillating conducting tip is used. 
Approach of the tip to the sample slightly changes 
the resonance frequency of the sensor, /o, and this 
frequency shift value, A/, can be used as a signal for 
the feedback loop to map the sample topography. The 
total interaction energy between tip and surface is a 
sum of the long-range electrostatic and van der Waals 
contributions and the short-range chemical interaction 
which provides the atomic contrast during scanning. 
As an add on to AFM the method of the Kelvin- 


probe force microscopy (KPFM) was developed 29 


that allows to obtain information about the local 
distribution of the contact potential difference between 
tip and surface ( Ulcpd )■ In this method the DC and 
AC (with the frequency w) voltages ( Udc and Uac, 
respectively) are applied between tip and sample and 
the measured first harmonic of the A/ signal is used 
to nullify the difference between Udc and Ulcpd- 
STM and AFM can be used for collecting the 
3D Ip{x,y,z ) and A f(x,y,z) sets of data. In the 
first approach the single ir(~) and A f(z) curves 
are collected on the dense grid and then they are 
combined in the 3D sets. In the second case the 
constant height 2D maps Ip{x,y) and A f(x,y) maps 
acquired with the usual scanning resolution at different 
equidistant distances and then combined in a 3D set. 
Such experiments are time consuming; therefore the 
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Figure 3. Methods of the graphene synthesis via (a) carbon segregation from the metal bulk, (c) CVD from hydrocarbons, and (c) 
carbon deposition from solid-carbon evaporation sources. 


careful postprocessing of the data, taking into account 
the thermal drift of the oscillating sensor and the 
scanner, is necessary. The detailed discussion of both 


approaches can be found in the literature 30-32 


2.3. NEXAFS and XMCD 

The method of near-edge x-ray absorption spec¬ 
troscopy (NEXAFS) and its respective extension, x-ray 
magnetic circular dichroism (XMCD), are used to ob¬ 
tain information about energy distribution of the un¬ 
occupied valence band states above Ep and for the 
determination of the magnetic moment of elements in 
the system [Fig. He)] [33, 
photon energy is scanned around the particular x-ray 
absorption threshold and the corresponding total elec¬ 
tron yield (TEY mode of NEXAFS) is measured as a 
drain current from the sample giving information about 
the absorption coefficient, which is proportional to the 
density of unoccupied states. If the photoelectron yield 
is measured by the electron multiplier then the low 
energy electrons can be removed from the signal by 
applying the negative potential of several volts to the 
electrostatic grid placed in front of the detector that 
allows to increase drastically the surface sensitivity of 
the method (partial electron yield or PEY). If linear 
polarized light is used in such experiments, then vary¬ 
ing the angle of the impugning light on the sample one 
can obtain the information about orientation of the va¬ 
lence band orbitals of the species adsorbed on the sur¬ 
face (the so-called search-light- like effect) [35]. If mag¬ 
netic sample is studied and circularly polarised light is 
used, then the absorption coefficient for the x-ray light 
depends on the relative orientation of the sample mag¬ 
netisation and projection of the spin of light on this di¬ 
rection (XMCD experiment). If two NEXAFS spectra 
for opposite directions of magnetisation are measured, 
then they can be used for the calculation of the spin 
(/is) and orbital (hl) magnetic moments of elements 


34 . In these methods the 


in the system |36,37|. 


2.J. Photoelectron spectroscopy 

In the method of photoelectron spectroscopy electrons 
are excited from the core levels or occupied valence 
band states and are analysed by their kinetic energy 
( Ekin) giving a replica of the electronic structure of 
solids [38]. Analysis of the core levels (x-ray photoelec¬ 
tron spectroscopy, XPS) gives an information about 
the chemical states of elements that might be influ¬ 
enced by the corresponding atomic coordination. In 
case of the angle-resolved photoelectron spectroscopy 
(ARPES) [Fig. He)] the electrons emitted from the 
valence band are additionally analysed by their emis¬ 
sion angle ( 9 ) and detected by the 2D detector giving 
photoemission intensity map for one fc-direction in the 
Brillouin zone. In the modern ARPES experiments the 
sample is usually placed on the goniometer having 5 or 
6 degrees of freedom. As shown in Fig.H c ); additional 
stepwise rotation of a sample around x-axis by angle 
/3 produces a stack of data I(Ekin,Q,0), which can be 
transferred to the respective set I(EB,k x ,k y ), where 
Eg and k x y are the binding energy and the respective 
component of the fc-vector of electron in the valence 


band 38-40. 


3. Graphene growth on metals 

There are three main methods of the graphene 
synthesis which are widely used now (Fig. [3j|: (a) 

segregation of the dissolved carbon from the bulk 
metal, (b) chemical vapour deposition (CVD) from 
hydrocarbons, and (c) deposition from solid carbon 
source (or molecular beam epitaxy, MBE). Here these 
methods are discussed in details. 
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Figure 4. (A) LEEM images of the graphene segregation on 
Ni(lll) at different stages: (a,b) the first layer segregates at 
1125K, (d-f) the second layer at 1050 K and (g-h) the third 
layer at 1050 K. Image (c) is a /rLEED pattern of the single¬ 
layer graphene/Ni(lll) demonstrating (1 X 1) structure. Data 
are taken from 4lJ with permission. (B) (a,b) Energy profile 
and the corresponding dynamic processes involved in graphene 
segregation, (c-e) STM images of graphene at different steps of 
segregation. Data are taken from [42] with permission. (C) (a) 
Real-time LEEM images of graphene segregation on Ru(0001) 
at 850° C. (b) Schematic illustration of the “downhill” graphene 
growth. Data are reproduced from |43| with permission. 


3.1. Segregation 


This method can be used for the graphene synthesis 
on metals which have high carbon solubility, like Ni 
or Ru 41 46 , but can work also for the metals which 


are heated to the high temperature that their carbon 
solubility becomes significant 47 51 . In this case, 


initially the metal bulk is loaded with carbon atoms 
at high temperature and then during cooling carbon 
atoms travel to the surface and form ordered layer of 
graphene. This method overcomes the main limitations 
of the CVD method as the multilayer graphene can be 
easily grown; thickness of a graphene layer depends 
on the cooling rate as well as on the loading and 
segregation temperature (Tg). 

For Ni the loading of bulk with carbon atoms is 


usually performed at temperatures above 1170 K 41 


Cooling of the loaded sample leads to the segregation 
of the carbon atoms at the surface and, depending 
on the segregation temperature and cooling rate, 
graphene of different thickness is grown. For single 
crystalline Ni(lll) the growth of single- and double¬ 
layer graphene was observed by means of AFM and 
Raman spectroscopy [52} . Similar studies performed 
by means of LEEM [Fig. [4j[ A)] shows that the first 
graphene monolayer is formed at Tg = 1125 K and 
the second layer grow at Tg = 1050 K 41 . It was 


also found that the third layer starts to segregate 
before second layer is formed. Compared to the 
single crystalline Ni surface, the segregation of carbon 
on polycrystalline films leads also to the growth of 
multilayer graphene which is attributed to the presence 
of grain boundaries in Ni that can serve as nucleation 
sites for multilayer growth [52] . 

The detailed mechanism of the graphene segre¬ 
gation on metallic surfaces (Ni, as an example) was 
presented in Ref. 42 and it consists of the follow¬ 


ing steps: (i) dynamic diffusion of the bulk dissolved 
carbon atoms with moving some of them to the sur¬ 
face, (ii) trapping of the surface-diffusing carbon atoms 
by defects or step edges on metal surface and cre¬ 
ation of the nuclei for the graphene growth, (iii) fur¬ 
ther graphene nucleation around the graphene centres 
[Fig- gB, a-b)]. The initial result of such process is 
shown in Fig. |4][B,c), where the start of nucleation of 
graphene is shown at several places. Decrease of the 
temperature leads to the rapid growth of graphene is¬ 
lands which finally form the complete graphene layers, 
that, as seen from the potential energy plot, is ener¬ 
getically favourable [Fig. |4][B,a)]. In such segregation 
process graphene is formed across the terraces, which 
are always present on surfaces of metals [Fig. |4](B,d-e)] . 

The validity of the above presented mechanism of 
segregation was also demonstrated for graphene growth 
on single-crystalline Ru(0001) [43,53 . Here initially 


the carbon atoms were absorbed into the bulk Ru 
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at 1150° C and cooling of the sample to 825° C leads 
to the segregation of the graphene on the Ru(0001) 
surface [Fig. QC)]. It was found that graphene islands 
grow via attaching of carbon atoms at the edges and, 
interestingly, it was found that graphene domains grow 
parallel to the substrate steps and across steps in the 
“downhill” direction. 


3.2. CVD 


Decomposition of hydrocarbons at high temperatures 
at the surface of metals is the easiest way to prepare 
continues graphene layers. This method applied to the 
polycrystalline metals, like Ni or Cu, demonstrated 
high perspectives to be used in graphene technology 
because huge high quality graphene layers of different 
thicknesses with the size of up to 30 inches were 
successfully synthesised (4[ [54j [55]. Several factors 
determine the growth of graphene on polycrystalline 
surfaces as pretreatment of the surface, density of the 
nucleation centers, flow of the hydrocarbon gases, use 
of additional flow gases (Ar, H 2 ), etc. 

Fig. 0 (A, upper row) shows secondary electron 
microscopy (SEM) images of graphene grown on 
unoxidized Cu foil at 1040° C for 0.5 hour and different 
flow rate of CH 4 [56]. The lower row demonstrate 
the dependence of the nucleation density as well 
as size of the graphene islands as a function of 
flow rate showing the optimal conditions for the 
graphene synthesis where nucleation rate and size of 
the graphene islands are balanced. Earlier it was shown 
that for the graphene grown on polycrystalline Cu 
two predominant planes were observed after graphene 
growth - Cu(100) and Cu(311) [57]. In this case 
graphene overgrow on these surfaces and facets and 
forms the complete carpet-like layer [Fig. 0 (B,a-c)]. 
However, the further experiments with the oxidation 
and subsequent reduction of Cu foils shows that 
predominant Cu(100) orientation might be obtained 
that can be used for the synthesis of the high-quality 
graphene layers [56 . 

The shape and size of graphene islands, which 
later coalescence in the complete layer, also depend 
on the ratio of H 2 or Ar to CH 4 58 59 . As was 


found in the recent systematic study, the higher ratios 
of Ar to H 2 result in graphene grains with a symmetric 
dendritic structures, whereas if ratio of Ar to H 2 is 
decreased, the more compact structures of hexagonal 


or circular shapes were observed 59 


The graphene growth on single-crystalline metallic 
surfaces is performed in UHV conditions that requires 
their careful preparation via cycles of ion-sputtering 
(Ne + , Ar + ) and annealing at high temperatures 60 


|67|. In several graphene/metal studies, when carbon 
has high solubility in the corresponding metal (Ni, Ru, 
or Rh), thick highly ordered films which were obtained 



Figure 5. (A) Upper row: SEM images of graphene islands 

on polycrystalline Cu. Numbers denote the flow rate of CH 4 . 
Lower row: the corresponding nucleation density and island 
size as a function of the CH 4 flow rate. Data are reproduced 
from [56] with permission. (B) STM images of (a) a large area 
of graphene on polycrystalline Cu, (b) flat area marked in (a) 
demonstrating atomic resolution of a graphene layer, and (c) a 
continuous graphene layer covering a monoatomic step of Cu. (d) 
Schematic presentation of graphene growth on Cu polycrystalline 
surface. Data are reproduced from [57] with permission. 


by e-beam evaporation of the respective material were 


used 68 77 


The close-packed surfaces Ni(lll) and Co(0001) 
have lattice constants which are very close to that 
of graphene (difference is of the order of 1 %) which 
make them very popular for graphene preparation 
and studies. Fig. [ 6 ] (A, upper part) 71 shows a 
series of XPS C Is spectra of graphene synthesised 
from C 3 H 6 on Ni(lll) at the respective temperatures 
(marked in the figure) demonstrating the optimal 
heating conditions: T > 650° C. At these parameters 
the intensity of the peak corresponding to the C 3 H 6 
fragments (marked by the black triangle at binding 
energy of 283 eV) vanishes completely. These results 
were later confirmed in STM studies, which results are 
compiled in Fig. [ 6 ] (A, lower part) 164], showing the 
high quality graphene grown on Ni(lll) at 650° C. 

In Refs. 78, 79 the XPS peak at 283 eV was 
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assigned to the carbon surface atoms located at step 
edges or other defect sites of the Ni(lll) surface (this 
peak also appears if the hydrocarbon source is replaced 
the solid carbon source). 

Later XPS and LEED experiments performed 
under in operando conditions |79] show that if clean 
Ni(lll) surface is exposed to hydrocarbons below 
500° C, then it leads to the formation of Ni 2 C which 
later converts to monolayer graphene via an in-plane 
mechanism (80 . Above 500° C graphene grows via 
replacement of Ni surface atoms leading to embedded 
epitaxial and/or rotated graphene domains 179 . 

The studies of graphene growth on Ru(0001) 
demonstrate that faceting of the metallic surface 
occurs with formation of large terraces where graphene 
forms a perfectly ordered layer 82 83 . Graphene 


in this system overgrows the metal steps and the 
process of growth is governed by the big Ru-atoms 
transport and these atoms migrate from “etched” 


terraces underneath a graphene layer 83 


Surprisingly, the studies of the graphene CVD 
growth on Rh(lll) showed that the temperature 
window for this process is very narrow: graphene 
was only found on the surface in the range between 
1016 K and 1053 K 


62 66 84 85 


Below this 

temperature range a rhodium carbide is formed. These 
results were compared with those obtained during 
temperature programmed growth (TPG) procedure 
(C 2 H 4 is adsorbed at room temperature and then metal 
heated to high temperature) and it was found that 
in this case the temperature range for the graphene 
formation is much broader between 808 K and 


1053K 84 


Formation of the graphene/Ir(lll) system was 
intensively studied in a series of LEEM and STM 
experiments 160, 86 - 89] . Fig. (g](B) shows a series of 
STM images collected after exposure of Ir(lll) to C 2 H 4 
at 1120 K for different time. Graphene coverage grows 
linearly in the beginning, asymptotically approaching 
full layer for the longer time. Compared to 
graphene/Ru(0001), on Ir(lll) a graphene layer grows 
in both directions - “uphill” and “downhill” across the 
metallic steps with preferential nucleation around the 
step edges, but growth on the flat terraces was also 
detected. Contrary to many other transitional metals, 
the growth of well-ordered graphene layer on Ir(lll) 
can be performed in very wide temperature range of 
970 K - 1470 K 60 . 


Compared to d-metals with open valence band 
shells, the CVD growth of graphene on the close d-shell 
metals is quite difficult. For example CVD graphene- 
growth on Cu(lll) is possible only at temperatures 
very close to the melting temperature. Fig. ( 6 (C) shows 
the STM results of such studies 81 . Here a graphene 
growth on Cu(lll) was performed at very high partial 


pressure of C2H4 (10 -5 mbar and dosing via nozzle) 
and cycles of fast sample flash-annealing to 1000° C, 
which is very close to the melting point of Cu (1083° C). 
In these experiments graphene layers of 0.35 — 0.8 
monolayer were studied. Later, the problem of the 
usage of the high temperature for graphene synthesis, 
which is very close to the melting point, was solved 
via irradiation of Cu(lll) and Au(lll) with ethylene 
(500 eV at 800° C) followed by the additional annealing 
at 900 - 950° C 90 . 


3.3. MBE growth 

This method is usually used for the preparation of 
the epitaxial graphene layers on the insulating and 


semiconducting surfaces at high temperatures 91-93 


as it was suggested that it can be easily adopted for the 
mass production and further graphene technologies. 
Recently it was shown that MBE growth can be 
successfully used for the preparation of graphene 
on Au(lll) in UHV conditions 94, 95] [Fig. [7](A)]. 
Carbon was deposited from an electron-beam-heated 
graphite rod at substrate temperatures ranging from 
770° C to 940° C. The quality of graphene/Au(lll) 
was monitored in situ by means of LEEM//xLEED 
and STM. On the initial stage of growth, graphene 
islands have a dendritic structure and, interestingly, 
in STM it was possible to resolve the untouched 
herringbone reconstruction of Au(lll) that indicates 
extremely small attraction energy between graphene 
and Au (< 13 meV/C-atom) [94]. MBE growth of 
graphene of different thickness (up to 3.5 monolayers) 
was also demonstrated on Ni(lll)/MgO(lll) [96| . 

Deposited carbon-contained molecules on the 
surface of metals might be considered as a solid source 


of carbon atoms for the synthesis of graphene [65 97 
For example, it was shown that Ceo can be used for 
the preparation of well-ordered graphene quantum dots 
of limited size on Ru(0001) using the high catalytic 
activity of this surface (97 [Fig. (7](B)]. If two precursors 
for the graphene growth are compared, C 2 H 4 and Cgo, 
then it was found that for the former one graphene 
is predominantly formed at the steps of substrate, 
whereas for the later case graphene dots have very low 
mobility at high temperature and might be randomly 
distributed over the Ru surface. 


4. Crystallographic structure of graphene on 
metals 

Synthesis of high-quality graphene layers on metals 
requires the perfect experimental conditions and 
always involves a high temperature. Any surface 
imperfections of metal or slight change of the synthesis 
conditions might lead to the appearence of different 
defects in graphene (vacancies, defect lines, etc.) as 
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Figure 6 . (A) Upper panel: (a) C Is XPS spectra of graphene grown at different temperatures on Ni(lll), (b) ratio of intensities 
of two peaks corresponding to C 3 H 6 fragments (F) and main graphene peak (G), (c) Doniach-Sunjic lineshape analysis of the C Is 
peak of the best-quality graphene. Data are taken from 71] w ith permission. Lower panel: STM images of graphene synthesised at 
different temperatures on Ni(lll). Data are taken from [64] with permission. (B) STM images of graphene layers synthesised on 
Ir(lll) after decomposition of C 2 H 4 (partial pressure 5 X 10 —10 mbar) at 1120 K for (a) 20 s, (b) 40s, (c) 160 s, and (d) 320s. Data 
are taken from [60] with permission. (C) STM results for graphene/Cu(lll): (a-c) large scale and (d) atomically-resolved small 
scale images. Data are taken from [81] with permission. 


well as to the different alignments of a graphene lattice 
on the surface of metal. Here we consider situations 
of the ideal alignment of graphene and metal close- 
packed surfaces when lattice vectors of both sublattices 
are parallel to each other, the so-called R0° structures 
(Fig.§. The discussion of other possible geometrical 
structures of graphene on metals as well as their 
interpretations can be found elsewhere 32.81,86 98- 


104 


Generally there is a lattice mismatch between 
graphene and metallic surfaces which is ranged 
from 1.3% for Ni(lll) to, e. g., 12.8% for Pt(lll) 


[Fig. IP ,b)]. Therefore during formation of graphene 
on Ni(lll) or Co(0001), graphene adopts the lattice 
constant of the metallic substrate and the so-called 
(1 x 1) structures are formed [Fig. |][a-c)] (although, 
as stated earlier, the slight variation of the synthesis 
parameters might lead to the formation of rotated 


structures, which were observed in the experiment 46 


74 105 ). For 4 d and 5 d metals adsorption of 


graphene on their close-packed surfaces always leads 
to the formation of the moire structures, which 
orientation depends on the relative orientation of the 
lattice vectors of graphene and Metal(lll) surface, 
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Figure 7. (A) (a) LEEM image (field of view: 5 /im) of 

graphene/Au(lll) grown by MBE at 950° C. The bright features 
are dendritic islands of graphene. The grey background is the 
bare Au surface, (b) LEED (40 eV) from an area of 2 fim in 
diameter, (c) STM images of graphene moire modified by the 
Au herringbone. Data are taken from [94] with permission. (B) 
STM images of graphene quantum dots on Ru(0001) synthesised 
from the pre-deposited C 60 molecules. Data are taken from i97 
with permission. 


whereas the corrugation of such structure is defined 
by the interplay between geometry and the local 
graphene/metal interaction around the high symmetry 
positions of the moire structure. 

Historically, the first method which gave informa¬ 
tion about crystallographic structure was LEED 12 - 


15 . Presently this method is used for the qualitative 


characterisation as well as for quantitative I — V LEED 


analysis of the graphene/metal structures 106-108 


However, the main method which provides an accurate 
analysis of the crystallographic structure on the large 
and small scale is STM, giving also access to the local 
electronic structure. 


(a) Top view 



(b) Side view 

ATOP HCP BRIDGE FCC 

• • •• •• •• • • •• • • • • .« • • •• •• grapheoe 


Metal (S) 
Metal (S-1) 


DOOOOOOOOOO Metal (S-2) 
(c) (1*1) structures 



Figure 8. Top (a) and side (b) views of a graphene layer on 
the close-packed (111) metal surface: (n X n)graphene/(m X 
m)Metal(lll) (here: n = 10, m. = 9). High symmetry 
adsorptions sites of the Metal(lll) surfaces are marked: 
circle (ATOP), triangle (HCP), rhombus (FCC), rectangle 
(BRIDGE), (c) Local high-symmetry (1 X 1) structures of the 
graphene/Metal(lll) interface. 


Fig. [9] shows a representative STM images of 
(a) (1 x l)gr/(l x l)Ni(lll), (b) (12 x 12)gr/(ll x 
ll)Rh(lll), (c) (13 x 13)gr/(12 x 12)Ru(0001), and 
(d) (10 x 10)gr/(9 x 9)Ir(lll). 

Graphene on Ni(lll) has a (1 x 1) structural 



Figure 9. STM images of (a) gr/Ni(lll), (b) gr/Rh(lll), (c) 
gr/Ru(0001), and (d) gr/Ir(lll). 
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A gr/Ru(0001) 


B gr/lf(111) 



Figure 10. (A) STM images of gr/Ru(0001): (a) and (b) are experimental and simulated STM images, respectively, obtained for 
the bias voltage Ut = -IV; (c,d) Simulated STM images for Ut = +1V and Ut = +3.5V. Insets show the apparent corrugation of 
graphene as a function of the bias voltage (experimental and theoretical data). Data are taken from [lQ9||llO| with permission. (B) 
Combined STM/AFM results for gr/Ir(lll): (a) and (b) are experimental and simulated STM images, respectively, where lower part 
of every image corresponds to Ut = —0.6 V and upper one to Ut = —1.8 V. (c) STM (top) and AFM (bottom) of gr/Ir(lll) where 
scanning mode was changed on-the-fly in the middle of the scan. Scanning parameters: Ut = 0.03 V, It = 1 nA, A/ = —475 mHz. 
(d) Extracted height profiles from STM and AFM data (see a-c). 


periodicity as the lattice constants of graphene and 
Ni(lll) are very close to each other [Fig. §a)]. 
Structural measurements based on the I — V LEED 
and photoelectron diffraction (PED) analysis give HCP 
(top — fee ) structure as the most stable arrangement 
of carbon atoms on Ni(lll) [Fig. | 8 fc)], with the 
distance between graphene and top Ni layer of 
2.11 A and very small corrugation of graphene of 
0.05 A 


106 


108 


These results were analysed in 
many DFT calculations (LDA, PBE, PBE-D, PBE- 
D2) and very good agreement was found between 
experiment and theoretical data 111-118 . The 


state-of-the-art calculations at the PBE-D2 level give 
the maximal bonding energy between graphene and 
Ni(lll) of about 160meV/C-atom 113 115 , which 


is much smaller than the lower limit for the chemical 
bonds, placing this value in the range of the van der 
Waals (vdW) interactions. In this system, graphene 
is strongly n-doped via electron transfer from the 
Ni4s on the graphene tt* states, that leads to the 
space-, energy-, and wave-vector-overlap of Ni3d and 
graphene 7 T states. As a result, several so-called 
interface states are formed, which are composed from 
Ni 3d XZ VZtZ 2 and C p z orbitals (see also discussion 
below) U][119|. Similar situation is also realised for 
the grapliene/Co(0001) interface 120 121 


Cu, 4d or 5d metals, always leads to the formation 
of the moire structures [see Figs. §C), 0 i 
Below we consider two interesting cases of graphene on 
Ru(0001) and Ir(lll), where correspondingly strongly 
and weakly corrugated graphene layers are formed. 

Graphene/Ru(0001) was investigated by differ¬ 
ent experimental methods, which give structural in¬ 
formation (corrugation of a graphene layer): I — V 
LEED (1 .5 A) 

(0.82 A) 


122 


123 


107 , surface x-ray diffraction (SXRD) 
, STM (between 0.5 A and 1.1 A de- 


Adsorption of graphene on other metals, like 


pending on the tunneling conditions) [109. 110||124- 

B- 

Presently, the widely accepted structure of 
graphene/Ru(0001) is (25 x 25)gr/(23 x 23)Ru(0001), 
meaning that 25 unit cells of graphene are aligned along 
23 unit cells of the Ru(0001) surface [107| 122 123] 
(see Refs. [129^1351 for the discussion of other smaller 
size models). These results were obtained in I — V 
LEED and SXRD experiments and compared with 
the DFT calculations within the GGA approach (PBE 
functional without inclusion of vdW interactions) 107. 
This modelling gives a corrugation of graphene of 
1.59 A (vs. 1.53 ± 0.2 A obtained in experiment) and 
correctly reproduce the obtained experimental data 
(I — V LEED curves, corrugation of the graphene 
layer, variation of bonds in the graphene layer, and 
corrugation of the topmost Ru layer). 
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However, the discussed model is quite expensive 
with respect to the computational time, and another 
model, namely (11 x ll)gr/(lOx lO)Ru(OOOl), was used 
for the investigation of the effect of vdW interactions 
on the electronic structure of gr/Ru(0001) 109 110 


by the parallel I — V LEED measurements 139 


In these works the PBE-D2 approach proposed by 
Grimme which includes vdW corrections 


are images as dark areas [Figs. |9[c) and 10 [B,a)] |67, 
138 . Increasing the bias voltage leads to the 
inversion of the imaging contrast and it becomes a 
direct one [Fig. [To[B,a)]. Similar to gr/Ru(0001), 
STM experiments on gr/Ir(lll) also show the bias 
dependence of the corrugation [Fig. [To|(B,d)], which 
can vary between 0.15 A and 1.0 A. Careful combined 
STM/AFM experiments of gr/Ir(lll) [32}|67 trace the 
changes in the imaging contrast and allow the accurate 
identification of the high-symmetry positions in this 


136 was 

used. This modelling shows that the corrugation of a 
graphene layer is 1.195 A (minimal distance be tween 
graphene and the top Ru layer is 2.195 A) [l09|, 
which is by 0.40 — 0.55 A smaller (depending on the 
structural model) compared to the values obtained 
without the vdW interaction 107| 129| 130| . Inclusion 
of the vdW interaction also leads to the dramatic 
increase of the bonding energy from 27meV/C-atom 
(no vdW) to 206 meV/C-atom (however as mentioned 
above this value is below the lower limit for the 
estimation of the covalent chemical bonding strength). 
The presented model correctly reproduces the STM 
data for gr/Ru(0001) [Fig. [l^A)] 109) 110 : STM 
topography as well as a value of corrugation (insets). 
For example, the inversion of the imaging contrast 
in STM was reproduced for the bias voltages around 
+2.5 V: above this point the topographically highest 
ATOP places of the gr/Ru(0001) structure are imaged 
as a dark areas in the STM data [Fig.[l0|(A,c-d)]. This 
effect was explained by the formation of the interface 
states at the HCP and FCC areas as a result of overlap 
of an unoccupied Ru(0001) surfaces resonance with the 
first image state component localized at these places of 
the moire structure [110+37 . 

Graphene on Ir(lll) forms a moire structure 
with a periodicity of (10 x 10)/(9 x 9) [Fig. |9jc)]. 
If low bias voltages between sample and tip are 
used for the imaging of this system (below ±1V, 
which is typical for STM experiments on metallic 
systems), then the inverted contrast in STM is 
observed, when topographically highest ATOP places 


Contrary to these results, the much larger corrugation 
of 0.6 A and 1.0 A for 0.39 ML and 0.63 ML graphene on 
Ir(lll), respectively, was obtained from x-ray standing 
wave (XSW) experiments [140| that was assigned to the 
increasing of the number of wrinkles on the surface. 

The graphene/Ir(lll) system was modelled within 
the framework of the DFT theory with vdW 
interactions accounted by means of a semiempirical 
DFT-D2 approach proposed by Grimme [136 or using 


a non-local correlation functional vdW-DF proposed 
by Dion et al. 141 143 These calculations lead to 


the following values for the graphene corrugation and 
the mean distance between the graphe ne l ayer and the 
top-most Ir layer: 0.31 A and 3.39 A 1 67 , 0.37 A and 
3.28 A 144], 0.35 A and 3.41 A [l40 . (The previous 
calculations performed at the GGA level without vdW 
interactions included give the minimal and maximal 
distances between graphene and Ir(lll) of 3.77A 
and 4.04 A, respectively [l45 .) All these simulations 
predict the relatively large distance between graphene 
and Ir(lll), which is very close to the distance between 
carbon layers in graphite. Analysis of the interaction 
shows that GGA approach leads to repulsion between 
graphene and Ir(lll) (a: +20meV/C-atom) and only 
inclusion of the vdW interaction gives the attraction 
in this system [140 j|. It was found that here graphene 
is slightly p-doped and there is a local overlap of 
the Ir5d.2 and C p z orbitals at HCP and FCC high- 
symmetry places of the graphene/Ir(lll) structure 
that leads to the appearence of the so-called hybrid 
states in the valence band of this system [67, 140 


Existence of such states was used for the explanation 
of the bias dependence of the change of the imaging 
contrast in STM data: ATOP positions are imaged 
as dark places at low bias voltages and become bright 
with increasing the voltage 67 . These results point 


the importance of the graphene-metal interaction in 
this system 67 101 , although in the beginning it 


was believed that graphene is nearly free-standing on 
Ir(lll) [146 . 

The relatively weak interaction between a 
graphene layer and Ir(lll) as well as between graphene 
and Pt(lll) leads to the fact that slight variation of the 
graphene synthesis parameters (temperature or/and 
partial pressure of hydrocarbons) might cause the ap¬ 
pearance of the rotational domains of graphene on the 
metallic surfaces, where graphene lattices are misori- 


structure [Fig. 10 B,c)]: ATOP positions are imaged 
as bright areas in the attractive regime (A/ < 0) of 
NC-AFM. The obtained corrugation of 0.3 A of the 
graphene layer in AFM experiment witht W-tip is 
very close to the theoretically calculated value (see 
discussion below). If CO-terminated tip is used, then 
corrugation of 0.42 — 0.56 A was obtained supported 


ented with respect to each other 32 81 86 98)104 


5. Electronic structure of graphene on metals 

Adsorption of a graphene layer on metallic surface 
always leads to the modification of the electronic 
structure of graphene. According to the modern state- 
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(a) gr/Ni(111), spin-up 


(c) gr/lr( 111) 




(d) gr/Cu(111); AE = 18 meV 



m •»— K —*• r 


(e) gr/lr(111); AE = 300 meV 



Figure 11. Calculated band structures of (a) graphene/Ni(lll) (spin-up) and (b) graphene/Cu(lll) along high-symmetry directions 
of the graphene Brillouine zone. Calculations are performed for (1 X 1) structures (small lattice mismatches are ignored in both 
cases). The it and a graphene-derived bands are highlighted by the red and blue colours, respectively, (c) Band structure of 
(10 X 10)graphene/(9 X 9)Ir(lll) unfolded for the graphene (1 X 1) primitive cell. Panels (d) and (e) show the zoomed regions of 
the respective electronic structures where energy gap for the graphene ir states is open. The corresponding energy gaps are marked 
in the figures. 


of-the-art electronic structure calculations graphene 
is always weakly bonded to metals: the bonding 
energy at the interface was found in the range of 
50 — 200meV/C-atom (see examples above), which is 
much lower than the lowest limit which is used for the 
description of chemical bonds. However, in all cases 
the original band structure of graphene around Ep 
(linear dispersion and the Dirac cone) is distorted via 
doping (n or p ) or/and via overlap with the valence 
band states of metal. The progress in the complete 
understanding of all changes in the electronic structure 
of graphene upon its adsorption on metals was reached 
in the recent publications 11 16 119 . The proposed 


universal model for the description of the graphene- 
metal interaction devides this process into two steps: 
(i) doping of graphene, n or p, by mobile s-electrons, 


that increases or decreases, respectively, the strength 
of the vdW interaction between graphene and metal 
and might lead, in the case of the n-doping, to the (ii) 
effective space-, energy-, and wave-vector-overlapping 
of the Cp z states of graphene and metallic d states. 
Depending on the relative energy-overlap of the n and 
d bands the Dirac cone in the electronic structure of 
graphene is fully destroyed (case of grapliene/Ni) or 
symmetry energy gap is opened directly at the Dirac 
point [Tb 119]. 

Fig. |11| shows calculated band structures of (a) 
graphene/Ni(lll) (spin-up) and (b,d) graphene/Cu(lll) 
for the HCP (top — fee ) arrangement (small lattice- 
mismatch is ignored) [Fig. | 8 ][c)]. The bonding ener¬ 
gies and the distances between graphene and metal 
surfaces are 67— 160 meV/C-atom and 2.11 A, respec- 
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a 


free-standing graphene 


a = 40° 



graphene/Ni(111) 


a = 40° 




Figure 12. A: comparison of the experimental and calculated NEXAFS spectra for (a) graphene and (b) graphene/Ni(lll) (a: 
is the incident angle of the linearly polarized light on the surface). Data are taken from 1 16| with permission. B: (a) CX-edge 
NEXAFS spectra, (b-e) LEED images, and (f) C Is XPS spectra of different graphene moire structures on 4 d and 5 d metals. Data 
are taken from [147] with permission. 


tively, for g r/Ni (lll) (depending on the used func¬ 
tional) 


113 


118 


148 


and 94 meV/C-atom and 3.02 A, 
respectively, for gr/Cu(lll) (PBE-D2). The result¬ 
ing band structure of these systems is very differ¬ 
ent. In case of graphene/Ni(lll), the initial n-doping 
of graphene shifts the Dirac cone below Ep, where 
energy- and wave-vector-overlap of the graphene tt 
states with Ni3d states, accompanied by the space 
overlap of the corresponding orbitals, leads to the for¬ 
mation of the so-called hybrid states at the interface 
between graphene and Ni(lll). The Dirac cone as well 
as original electronic structure of graphene in the vicin¬ 
ity of Ep are fully destroyed 116,117 . If graphene is 
adsorbed on Cu(lll), then in the first step the doping 
shifts the graphene Dirac cone below Ep, but over¬ 
lap of the graphene 7 T and Cu3d states appears at the 
large binding energies below the Dirac cone. The deep 
analysis of this system shows that the so-called symme¬ 
try energy gap is opened directly at the Dirac cone of 
graphene due to the overlap of the C p z states from 
two carbon sublattices of graphene with the Cu3d 
states of different symmetry from the top-most Cu 
layer [l 19 1491. 

Fig.|ll[c,e) shows calculated band structure of the 
(10 x 10)graphene/(9 x 9)Ir(lll) system unfolded to 
the graphene (1 x 1 ) primitive unit cell according to 
the procedure described in Refs. 1150, 1511 with the 


band structure have a spaghetti- style, one can clearly 
recognise the band dispersions of the graphene-derived 
n and er states. These calculations show that graphene 
in this system is p-doped and the Dirac point is located 
at 150 meV above Ep and the energy gap of ~ 300 meV 
is opened for the n states at the K point. This result 
is in good agreement with recent ARPES data for 
graphene/Ir(lll) 101 ], where lowering of the Dirac 


point was induced by the K adsorption that allows 
observation of the energy gap of ~ 100 meV at the 
A'-point. 

Core-level based electron spectroscopies (XPS 
and NEXAFS) are powerful methods for investigation 
of the electronic structure of the graphene-based 
systems. Fig. [12] shows the NEXAFS and XPS spectra 
of different graphene-metal systems in comparison 
with those for graphite (con sists of the layers of 
graphene separated by 3.4 A) |l!2 


113 

116 

147 

152 


153 . NEXAFS spectra of graphite are very good 


example for the demonstration of the search-light -like 
effect as the relative intensities of the Is —>• n* and 
Is —* a* absorption bands are changed if the incident 
angle a is varied 


152 


code BandUP Em Although, from the first look the 


NEXAFS and C Is spectra can be used for the 
investigation of the orbital overlap of the valence band 
states of graphene and metallic substrate. In case 
of graphene on Ni(lll), Rh(lll), and Ru(0001) the 
strong modification of the Is — > tt* absorption band is 
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Figure 13. ARPES intensity maps for (a) graphene/Ni(lll) and (b) graphene/Ir(lll) measured along the V — K direction of the 
Brillouine zone, (c) Zoom of the ARPES map of graphene/Ir(lll) marked by the dashed line in (b). Photon energy is hv = 65eV. 


observed (Fig. 12) compared to the one for graphite, 
indicating the strong gr7r - Metal d overlap at the 
interface with a formation of hybrid states, whereas 
NEXAFS spectra for gr/Pt(lll) and gr/Ir(lll) are 
moderately influenced. [The respective changes in the 
Is — a* absorption band for gr/Ni(lll), gr/Rh(lll), 
and gr/Ru(0001), compared to HOPG and gr/Pt(lll), 
are related to the partial sp 2 -to-sp 3 re-hybridization 
for the former systems 74 147 .] C Is XPS spectra 
for these systems are different reflecting the influence 
of the crystallographic structure of the system on 
the electronic structure. In case of gr/Ni(lll) this 
spectrum consists of one line (shifted by ~ 0.4 eV 
with respect to the spectrum for graphite due to the 
charge transfer) indicating that graphene on Ni(lll) 
is flat [see STM image in Fig. |9][a)] and that the 
inequivalency in position of two carbon atoms in 
the graphene lattice has extremely small influence 
on the electronic structure. C Is spectra for the 
graphene moire structures are different. Following the 
previous considerations, the mean distance between 
graphene and Ir(lll) or Pt(lll) is quite large of 
ss 3.4 A with a graphene corrugation of about 
0.3 A, and the C Is spectra for these systems have 
only one component. For the gr/Rh(lll) and 
gr/Ru(0001) systems, graphene is strongly corrugated 
(« 1.1 A) with places where the distance between 
graphene and metallic surface is close to the one 
for gr/Ni(lll) (« 2.1 A) and where strong orbital 
overlap occurs. Therefore the C Is XPS spectrum 
for these graphene-moire structures consists of two 


components corresponding to valleys (high binding 
energy component C2) and hills (low binding energy 
component Cl). 

Strong orbital intermixing of the graphene 7r and 
Ni3d valence band states at the interface, discussed 
above, leads to the appearence of the induced magnetic 
moment of carbon atoms which was detected by means 
of XMCD 112|152|[l54 . This observation is supported 
by the electronic structure calculations of the induced 
magnetic moments 111 112 as well as by the 


spin-resolved ARPES measurements of the exchange 
splitting of the 7r band in the graphene/Ni(lll) 
system |74j. 

Fig. |T3] shows ARPES intensity maps mea¬ 
sured for two systems, graphene/Ni(lll) (a) and 
graphene/Ir(lll) (b-c), which can be considered as two 
limit cases of the graphene-metal systems where strong 
and weak overlaps, respectively, of the valence band 
states at the interfaces are observed. 

Grapliene/Ni(lll) interface has (1 x 1 ) structure 
and strong orbital mixing of the graphene n and Ni 3 d 
valence band states is found at the interface 74 


11 l|[TT3 155] [Figs. [TT]a) and [E^a)]. In ARPES 


maps one can clearly identify graphene 7r and er 
states which are shifted to higher binding energies 
(compared to free-standing graphene) due to the 
electron doping. The effect of orbital mixing is 
clearly visible around the A'-point of the Brillouine 
zone where several hybrid states are formed in the 
energy range between E F and 2.3 eV. These two effects, 
doping and states-overlap, shift 7r states by ~ 2.4 eV 
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to higher binding energies, compared to « 1 eV for 
the <r states, which are influenced by the doping 
only. Valence band states of Ni are also strongly 
modified as compared to the clean Ni(lll) surface do 
to the charge redistribution at the interface |69 . The 
similar behaviour of the valence band states is also 
found in other ARPES studies of the graphene-metal 
systems, where significant overlap of the valence states 
occurs at the interface: graphene/Co(000 1) [156], 
graphene/Rh(lll) [157] , graphene/Ru(0001) [131| 158 
|159| , graphene/Re(0001) [160| . 

The recent low-temperature ARPES experiments 
(40K) on graphene/Ni(lll) and graphene/Co(0001) 
were focused on the investigation of the energy region 
where hybridization of Ni 3 d and graphene n states 
These results suggest the existence of 


occurs 


164 


intact Dirac cones in these graphene-metal systems via 
over-doping of graphene which is contact with Ni (Co). 
This consideration contradicts the present description 
of these systems and these observations require further 
experimental and theoretical analysis. 

The electronic structure of graphene/Ir(lll) is 
shown in Fig. 13 b-c), where one can clearly identify 
graphene-derived tt and a bands [10T| p6l[l63l [165] - 


167 . In this system graphene is p-doped and Dirac 
cone is found ss 100 meV above Ep. Due to the weak 
interaction between graphene and Ir(lll), the energy 
bands of graphene are not disturbed via intermixing 
with Ir bands. This effect is clearly indicated by 
the fact that the Rashba-split Ir p z surface resonance 
(SR) around the T-point remains intact upon covering 
Ir(lll) surface with graphene (only upward shift was 
detected) (166,[167| (confirmed by the XPS studies 
of the surface core-level shifts between Ir(lll) and 
gr/Ir(lll) [168] ). As was shown above, graphene 
on Ir(lll) forms a moire structure with periodicity 
(10 x 10)gr/(9 x 9)Ir(lll) [Figs. [8](a), ]9|d), and 10 


and effect of this additional periodicity with small 
modulating potential is visible in ARPES maps as 
additional photoemission replicas for the main 7r 
and a emission bands [Fig. 13[b-c)[. Intersections 
of these replica bands and main bands produces 
the so-called mini-gaps in the electronic structure of 
graphene due to the avoid-crossing mechanism. For 
example such replica bands and the corresponding 
mini-gaps at the binding energies of sa 1 eV and 
ss 2.5 eV are clearly visible in Fig. |T3[c). The 
similar behaviour of the valence band states of 
graphene was found also for other graphene-metal 
systems with weak overlap of the valence band states 
at the interface: graphene/Cu(lll) [l04|[l69l[T70| 


graphene/Pt(lll) [47 171| , graphene/Au(lll) [95|. 


6. Graphene hetero- and nano-structures on 
metals 

The electronic structure of graphene on metals can be 
tailored in different ways that might help to understand 
the effects observed in the real electron- and spin- 
transport graphene-based devices. There are many 
ways for such modifications and several of them will 
be considered here: (i) intercalation of different species 
between a graphene layer and metallic substrate, (ii) 
adsorption of atoms, molecules and clusters on top 
of the graphene-metal system, and (iii) preparation 
of graphene objects of reduced dimensionality (flakes, 
quantum dots, nanoribbons). 


6.1. Intercalation 


There are many recent experimental and theoretical 
works where different species were intercalated between 
graphene and metallic substrate with the aim to per¬ 
form controllable modifications of the crystallographic 
and electronic structures of the system. Among these 
materials are: (i) alkali and alkali-earth metals, which 
induce n-doping of graphene 155 ,162 172 175], (ii) 
sp-, d-, and /-metals, where doping depends on the 
intercalated metal [6^ [TI^[n^[TI^[T44l[T49l|157[|159 


160 174 176-204 


and (iii) molecular species, which 
decouples graphene from the metallic support, like oxy¬ 
gen [163[[205H209] , CO [210], and C 6 


^60 


2J f 213 . 


The intercalation of monovalent species (alkali 
atoms and halogens) as well as oxygen and CO in 
graphene/Ir(lll) was studied theoretically in Ref. [161 
[Fig. [3 A)]. It was shown that these materials 
effectively decouple a graphene layer from the metallic 
substrate and induce the corresponding doping of 
graphene. For example alkali metals induce n- 
doping and oxygen or halogens make graphene p- 


doped. Authors of Ref. 161 found a correlation (linear 


dependence) between bonding of graphene to substrate 
and the doping level of graphene [Fig. [SA, a -b)]. 
This fact indicates that for the studied species (simple 
case of the charge transfer) the interaction between 
graphene and intercalants has ionic character. Authors 
also found a linear dependence of the shift of the C Is 
line on the doping level. 

The results of this theoretical work are con¬ 
firmed by the experimental ARPES data shown in 


Fig. [M] for (B) graphene/Cs/Ir(lll) [162] and (C) 
graphene/0/Ir(lll) [163 . One can clearly see that 


in both cases graphene is effectively decoupled from 
Ir(lll) as the corresponding replica bands induced 
by the moire structure of gr/Ir(lll) are not vis¬ 
ible in the respective intercalation systems. The 
obtained doping level of E — E F = — leV for 
gr/Cs/Ir(lll) Pig- 0B)] [162| and E — Ep = 
+0.64 eV [Fig. [14JC)] [163] correspond to the strongly 
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(a) ,-. (b) (c) 




Figure 14. A: (a) Partial C p z density of states for different intercalation graphene/X/Ir(lll) systems as compared to the free¬ 
standing graphene, (b) and (c) Calculated correlations between bonding energy of graphene and the corresponding core-level shift 
of the C Is line, respectively, in the intercalation systems from (a) and the corresponding doping level. Data are taken from 161 
with permission. B: ARPES maps measured along the direction perpendicular to T K of the Brillouine zone of graphene for (a) 
graphene/Ir(lll), (b) graphene/0.5 ML Cs/Ir(lll), and graphene/1 ML Cs/Ir(lll). Data are taken from 1 62] with permission. C: 
ARPES maps (a-c) and the corresponding constant energy cuts measured for graphene/Ir(lll) before and after intercalation with 
oxygen. Data are taken from 163] with permission. 


n- and p-doped graphene, respectively, and these val¬ 
ues are in very good agreement with the doping levels 


obtained in DFT calculations [Fig. 14 [A,a)] 11611. 

Intercalation of graphene on the particular 
metallic substrate might help to synthesise the artificial 
system which cannot be prepared in a direct way. 
Such examples are shown in Fig. [l5|(A,B). In the 
first case graphene/Ir(lll) was intercalated with Ni 
and strongly corrugated graphene, compared to weakly 
buckled graphene/Ir(lll), was obtained 144 . Nickel 


layer was deposited at room temperature on top 
of graphene/Ir(lll) and then system was annealed 
at temperatures up to 800 K. It was found that 
Ni intercalates under graphene as shown by STM 
and core-level XPS [Fig. |l5|[A,a-c)]. In this case 
Ni grows in a pseudomorphic way at the interface 
between graphene and Ir(lll) that periodicity and 
symmetry of the system is preserved [LEED and 
STM in Fig. 15 [A,a-c)]. As was found in STM 
the imaging contrast in this case is changed and 
STM images are very similar to those for the 
strongly buckled graphene/Ru(0001) system. In the 


same experiments with identical imaging conditions 
(bias voltage and tunneling current), the corrugation 
of 0.6 A was measured for gr/Ni/Ir(lll) compared 
to 0.25A for gr/Ir(lll). DFT simulation of 
gr/Ni/Ir(lll) yields the corrugation of graphene of 
1.51 A [1.195 A for gr/Ru(0001)] with a minimal 
distance between graphene and Ni-layer of 1.94 A 
(2.195 A for gr/Ru(0001) and 2.1 A for gr/Ni(lll)). 
It is interesting that more than 70% of carbon atoms 
are placed at the distance 2.0 — 2.2 A above Ni layer 
giving as a results the STM images similar to those for 
gr/Ru(0001). 

XPS spectra measured before and after Ni- 
intercalation in gr/Ir(lll) demonstrate the strong 
modification of the C Is line [Fig. [l5](A,c)] : the single 
peak at 284.1 eV for gr/Ir(lll) is shifted to 284.9eV 
for gr/Ni/Ir(lll) and shows a strong asymmetry after 
intercalation indicating the existence of two regions in 
a graphene layer where weak and strong orbital mixing 
are observed [similar to gr/Ru(0001) where double 
peak structure of C Is was found, Fig. HlB,f)]. 

The band structure of the graphene/Ni/Ir(lll) 
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Figure 15. A: (a-b) Large and small scale STM images of graphene/Ni/Ir(lll) demonstrating the process of intercalation and the 
structure of the system. Inset of (a) shows the corresponding LEED. (c) Comparison of C Is XPS spectra at different stages of 
intercalation of Ni in graphene/Ir(lll). (d) ARPES map for graphene/Ni/Ir(lll) measured along the P — K direction. Data are taken 
from 144 with permission. B: Large and small scale STM images of graphene/Cu/Ir(lll). Inset of (a) shows the corresponding 
LEED. (c) Comparison and the respective fit of C Is and Ir4/ XPS spectra of graphene/Ir(lll) and graphene/Cu/Ir(lll). (d) 
ARPES map on the large energy scale and (e) the corresponding zoom to the energy range close to Ep for graphene/Cu/Ir(lll). 
Data are taken from 149] with permission. 


system measured by means of ARPES is shown in 
Fig. [T5j A,d) and it is very similar to the one for 
graphene/Ni(lll): the graphene-derived 7r states are 
found at ss 10 eV of the binding energy at the F-point 
and the Ni3d weakly dispersing bands are located just 
below Ep- The position of the n band at the Ap¬ 
point depends on the thickness of the intercalated Ni 
layer: it is changed from 2.16 eV for 1 ML-intercalated- 
Ni to 2.65 eV for multilayer-intercalated Ni 144] (this 
effect is connected with the narrowing of the d-band 
with decrease of the thickness of the Ni layer and 
demonstrate the effect of Ni3d - graphene n orbital 
overlap or hybridization). These results support the 
model discussed above (see Sec.[5]and Ref. 119 ) where 
modification of the electronic structure of graphene on 
the open d-shell metal depends on the relative energy 
position of the shifted Dirac cone and the d-band: 
complete space-, energy-, and wave-vector-overlap 
leads to the destroying of the electronic structure of 
graphene around Ep with formation of the so-called 
hybrid states below Ep. 

Another example of the artificial system, graphene/ 
Cu/Ir(lll), obtained after intercalation, is shown in 
Fig. [149 . Similar to the previous system, 

the copper layer grows pseudomorphically underneath 
graphene as deduced from LEED and STM data 


tained intercalation system gr/Cu/Ir(lll) was always 
imaged in STM in the true topographic contrast and 
its variation upon changes of the bias voltage was not 
detected. This system was modelled by DFT (PBE- 
D2) and these calculations gave very good agreement 
with the experimental STM images. Relaxation of the 
geometry of graphene/Cu/Ir(lll) yields a corrugation 
of graphene of 0.229 A with a mean distance between 
graphene and Cu layers of 3.02 A. These calculations 
also predict the ?r-doping of graphene with the position 
of the Dirac point at 0.45 eV below Ep. 

XPS data for gr/Cu/Ir(lll) confirm these pre¬ 
dictions of the DFT calculations and consistent with 
STM data [Fig. ]l5][B,c)]. The Cls XPS line is shifted 
to the larger binding energies and can be fitted with 
two components which ratio gives the ratio of the ar¬ 
eas for hills and valleys of the gr/Cu/Ir(lll) struc¬ 
ture. The respective changes are also observed for the 
Ir4/ XPS line: the emission from the interface com¬ 
ponent (i) is suppressed for gr/Cu/Ir(lll) compared 
to gr/Ir(lll) and its binding energy is increased indi¬ 
cating that chemical surrounding for these Ir atoms is 
also changed. 

ARPES intensity maps for graphene/Cu/Ir(lll) 


are shown in Fig. 15 B,d-e). They show an n-doping 


[Fig. 15 [B,a-b)]. Contrary to graphene/Ir(lll), the ob¬ 


of graphene in this system of about 0.688 eV. Also 
the clear hybridization between Cu 3d and graphene 7r 
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bands is visible in the energy range of ~ 2 — 4eV 
that leads to the opening of the energy gaps due to 
the avoided-crossing mechanism. Surprisingly, a large 
energy gap of 0.36 eV is opened directly at the position 
of the Dirac point. DFT calculations (PBE-D2) give 
a value of 0.15 eV. Similar behaviour of graphene 
valence band states and gap opening at the Dirac 
point was also observed in Refs. 160, 186 for other 
graphene-metal intercalation systems (Cu, Ag, Au as 
intercalants). The proposed in Ref. [186] explanation 
that value of the gap depends on the doping level of 
graphene does not supported by other ARPES data 
(for example, the absence of the energy gap for alkali- 
metal intercalated graphene-metal interface 162||l75 


hydrogen (these maps show the electronic structure of 
graphene measured along the direction perpendicular 


Authors of Refs. 119 149 propose a model that 


appearing of the energy gap at the Dirac point in the 
electronic structure of graphene, adsorbed on the open 
d-shcll metal, is due to the violation of the symmetry 
of the electronic states as the p z orbitals of carbon 
atoms from different sublattices of graphene overlap 
with the d orbitals of the different symmetries of the 
same metal atom below a graphene layer. This model 
was tested for different doping levels of graphene and 
different distances between graphene and surface of d - 
metal and clear correlation between these parameters 
and the width of the energy gap was found 119 . These 


will be modified 214 -221 


of adsorbate and valence states of graphene, and in 
ideal case the violation of the sublattice symmetry 
in graphene that can lead to the opening of the 
energy gap at the Dirac point. Secondly, if e. g. 
moire graphene structures on metals are studied, then 
deposition of molecules or metals can lead to the 
formation of ordered arrays of molecules or clusters 

In this case the 


61 134 145 165 222-244 


atomic-H on graphene/Ir(lll) [214 215 . The panels 


to T — K, see also Fig. 13 b,c)). The clear modification 


changes have a huge impact on the electronic structure 
of graphene around the Dirac point and might influence 
the transport properties of graphene at the graphene- 
metal interfaces. 

6.2. Adsorption on graphene/metal 

The aim of the studies of adsorption on top of 
graphene is manyfold. Firstly, one can expect the 
selective adsorption of atoms of different nature with 
the result that the electronic structure of graphene 

doping, overlap of orbitals 


of the emission picture is visible which indicated by 
opening of the energy gap at Ep for higher exposure 
time, broadening of the emission lines, and reduction 
of the photoemission intensity. The parallel STM 
measurements performed in the similar experimental 
conditions [Fig. [l6]( A,d-i)] demonstrate the selective 
hydrogen adsorption at the FCC high-symmetry places 
of the graphene/Ir(lll) moire structure on the initial 
stages of the experiment; increase of the adsorption 
time leads to the formation of the elongated structures 
as observed in STM. However, even long hydrogen 
dosing allows to identify that adsorption of H-atoms on 
graphene/Ir(lll) follows its moire structure [see FFT 
analysis in Fig. [l6|[A,i)]. 

DFT analysis of the hydrogen adsorption on the 
graphene/Ir(lll) structure performed in the same 
work E3 shows that it is energetically favourable 
for hydrogen to form a “graphane-like” islands at the 
FCC places of gr/Ir(lll), where one of the C-atoms 
from the unit cell of graphene (C -hep) is bonded to the 
H atom and the second C-atom (C -top) is bonded to 
the underlying Ir atom. Such adsorption configuration 
leads to the local rehybridization of the carbon atoms 
in the graphene layer from sp 2 to sp 3 . Such regions 
of the “graphane-like” structures have very large band 
gaps (might reach 3.5 eV for true graphane 245 ). In 
reality, DFT calculations show that, e. g. at 23% 
coverage of graphene with atomic hydrogen, a band 
gap of 0.73 eV can be obtained. As was concluded on 
the basis of the experimental and theoretical data, the 
band gap opening in the H-graphene/Ir(lll) system is 
caused by the confinement effect in the residual bare 
graphene regions and the broadening of the bands is 
due to the increased uncertainty in the wave-vector for 
electrons in such structure 


214 


The investigations of the cluster formations on top 
of graphene/metal moire structures and their influence 
on the electronic structure of graphene were performed 
in a series of experimental and theoretical works listed 
above. The first experiments were carried out for 
Ir-clusters on graphene/Ir(lll) [Fig. 16 [B,a-c)] 145 


on top 

properties of the complete system can be modelled on 
the basis of its single element - single cluster - and 
different properties of such objects (catalytic activity 
or magnetic properties for magnetic elements) can be 
studied in details. 

Figure [l6|[A) shows an example of adsorption of 


(A,a-c) show the evolution of the electronic structure of 
graphene around the AT-point of the Brillouine zone as 
a function of the exposure time of the system to atomic 


These STM results combined with DFT calculations 
show that Ir atoms preferentially nucleate at the 
HCP high-symmetry regions of graphene/Ir(lll) 
forming a regular arrays of Ir clusters. According 
to DFT calculations adsorption of Ir atoms on 
graphene/Ir(lll) leads to the local rehybridization 
(from sp 2 to sp 3 ) of carbon atoms at the FCC and 
HCP regions (similar to discussed earlier). In this 
case the C-atom which is placed at the hep or fee 
position is covalently bound to the Ir atom placed 
above the graphene layer and the C-atop at the top 
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Figure 16. A: (a-c) ARPES maps recorded along the A — K — A direction (see inset of (a)) for Osec, 30sec, and 50sec, respectively, of 
dosing of atomic hydrogen on graphene/Ir(lll). (d-h) STM images of graphene/Ir(lll) taken after very small time, 15sec, 30 sec, and 
50 sec, respectively, of dosing of atomic hydrogen, (i) FFT analysis of image (h). Data are taken from 2 |T| with permission. B: STM 
images of (a) graphene/Ir(lll) and (b) 0.02 MLIr/graphene/Ir(lll). (c) Structural model of graphene/Ir(lll) with C-Ir distances 
obtained from GGA calculations, (d) Constant energy ARPES maps for graphene/Ir(lll) and 0.15 MLIr/graphene/Ir(lll). (e) 
Scheme of the Dirac cones for the main photoemission band and its replicas, (f) Colour coded scheme of the distribution of the 
moire potential in the systems. Data are taken from |l45| 165] with permissions. 


position is covalently bound to the underlying Ir 
The similar situation is observed 


145 222,226 


atom 

for other metallic adsorbates 224 231| 243 

Adsorption of clusters or molecules on the 
graphene/metal systems can lead to the changes in 
the electronic structure of a graphene layer (change 
of the doping level or/and modification of the energy 
dispersion of the graphene valence band states) as 
was demonstrated in the experiment. The electronic 
structure of the Ir/graphene/Ir(lll) system was 
studied by means of ARPES and these results are 
shown in Fig. [l6](B,d-f) 165|. As was found in the 
experiment the deposition of Ir on gr/Ir(lll) (the 
nominal thickness of Ir was 0.15 ML that corresponds 
to one 13-atom Ir cluster per HCP region of gr/Ir(lll)) 
leads to the opening of the energy gap in the electronic 
structure of graphene at the .A-point and position 
of the lower part of the Dirac cone is placed by 
200 ± 20meV below Ep. The important result is that 
the symmetry of the system is reduced from six-fold to 


three-fold as can be deduced from the comparison of 
photoemission maps for gr/Ir(lll) and Ir/gr/Ir(lll) 
[Fig. 16[B,d-e)]. This reflects the fact that local 
rehybridization from sp 2 to sp 3 appears at the HCP 
regions upon Ir clusters formation and this perturbs 
the photoemission intensity suppressing three out of 
the six photoemission replicas |165| . The formation of 
the Ir arrays of clusters on gr/Ir(lll) also increases 
the width of the mini-gaps where replica bands cross 
the main graphene 7r-band from 240 ± 20meV for 
gr/Ir(lll) to 330±20meV for Ir/gr/Ir(lll) indicating 
the effect of the strengthening of the modulating moire 
potential. As a result the strong anisotropy of the 
group velocity of the graphene 7r states along V — K 
and perpendicular to this direction was measured that 
indicates the importance of the periodic potential on 
the transport properties of graphene and can help 
to tailor the transport properties of graphene in the 
future devices. The discussed effects (doping and 
increasing of the width of mini-gaps) were dramatically 
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increased by the co-adsorption of Ir and Na atoms 
on graphene/Ir(lll) [ 233] , that was explained by the 
even stronger modulations of the moire potential in the 
obtained systems. 

Magnetic properties of magnetic clusters on top 
of graphene/Ir(lll) were studied by means XMCD 
in Ref. [228| . Different clusters were studied: 
Pf 13 C 026 , Pti3Fe26, Iri3Co26, Ir5oC°500) Co27oo- The 
morphology and quality of the studied systems were 
verified by STM. It was found that for small clusters, 
e. g. Pti 3 Co 26 , magnetization curves do not reach 
the saturation even at magnetic fields of 5 T and no 
hysteresis down to 10 K was observed. The extracted 
spin and orbital magnetic moments of Co in such 
clusters are p, s = 1.5 ± 0.2 p,B and /q = 0.22 ± 0.02/is, 
respectively (fj, s = 1.62/iq and /q = 0.15 /j,b for 
bulk Co). Larger clusters (Ir 5 oCo 5 oo and C 02700 ) 
demonstrate the non-zero coercivity, which vanishes 
around 40 K. As was found, the spin (orbital) magnetic 
moments for such clusters are slightly increased 
(decreased) by 0.2/iq (0.02 — 0.04 pb) compared to the 
values for small clusters. 

In Ref. [246] the regular array of Pt clusters was 
formed on graphene/Ir(lll). This system was exposed 
to CO gas and its stability was tested by means 
of STM, XPS, and results were analysed with DFT 
calculations. It was found that for the clusters of the 
size of few tens of atoms such adsorption of CO leads to 
the sintering of the clusters via Smoluchowski ripening 
- cluster diffusion and coalescence. Larger clusters 
upon exposure to CO remain stable but form three- 
dimensional larger agglomerates. These effects were 
explained by the weakening of the Pt-C interaction 
upon CO adsorption in the 2-fold edge bridge positions 
of the cluster between Pt atoms. As was found such 
position is energetically more favourable compared to 
the 1-fold position when CO is bonded to every Pt 
edge atom of the cluster. Such 2-fold adsorption of 
CO leads to the increase of the Pt-C distance by 0.9 A 
and clusters become less bonded to the gr/Ir(lll) 
substrate that increase their mobility and probability 
to coalescence to large structures. Such studies of 
the stability of the cluster arrays on graphene/metal 
systems might shed light on the understanding of their 
catalytic properties in future studies. 

6.3. Graphene nanoohjects: nanoribbons, nanoflakes, 
quantum dots 

A graphene nanoribbon (GNR) is a narrow strip 
of graphene, which structure and the electronic 
properties are defined by the edge morphology and 
the width. The morphology is defined by the chiral 
indexes (n, m) (similar to CNTs) or chiral angle 
6. GNRs can be classified, with respect to the 
morphology of edges, as armchair, zigzag or chiral. 



Figure 17. A: STM images of (a) 7-AGNRs and (b) chevron- 
type GNRs on Au(lll). In both panels the STM images are 
overlaid with DFT-based simulations of STM (grey) and atomic 
models of GNRs (C: blue, H: white), (c) 3D representation 
of STM image of 7-AGNR. (d) STS spectra of 7-AGNR (red) 
and clean Au(lll) (black), (e) STM profile across 7-AGNR. 
(f) Series of colour-coded STS spectra taken across 7-AGNR. 
Data are taken from [247[|248| with permissions. B: (a) STM 
image of unidirectionally aligned 7-AGNRs on Au(788). (b) 

Structure of hydrogen-terminated 7-AGNR and relevant lattice 
parameters, (c) Brillouin zone of graphene (black hexagon) 
and the one- dimensional Brillouin zone of 7-AGNR (red), (d) 
ARPES intensity plot I(E — Ep,k ||) recorded along the ribbon 
axis. Data are taken from [248] with permissions. 


The electronic structure of GNRs is different from the 
one of graphene due to the confinement effects and 
as was shown theoretically, depending on the edge 
morphology and the width, GNRs can be metallic, 
semimetallic or semiconducting 249 . According to 


these considerations 249-251 , zigzag GNRs (ZGNRs) 
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show a sharp peak in the density of the electronic states 
at E F that leads to a net spin polarization at the edges; 
antiferromagnetic coupling between two edges opens 
a small fundamental gap. The armchair graphene 
nanoribbons (AGNRs) can be either semimetallic 
or semiconducting depending on the width of the 
nanoribbon [252 256 . 

However, the contact of graphene nanoribbon with 
a metallic surface can drastically change the electronic 
structure of GNR due to the effects of doping and 
orbital overlap of the electronic states of nanoribbon 
and metal as was shown above for a graphene-metal 
interfaces. Thus, the previous theoretical consideration 
have to be revisited with the aim to account for 
the metallic contact to GNRs. For example in 
Ref. 257 it was shown that although the interaction of 


247 


precursors 

these nanoribbons give a bandgap of 2.3 ± 0.1 eV for 


7-AGNR [see Fig. [L7jA,c-d)] [248] and 1.4±0.1eV for 
13-AGNRs 268 . These values can be compared with 


7-AGNRs and chevron-type GNRs, respectively 270 


from GNRs on Au(lll) that leads to the surface po¬ 
larization, that influences the width of the energy gap 
of GNRs. 

Later the above presented method of synthesis of 
GNRs 247 was used in Ref. 267 where structural 


nanoribbons with noble-metals is weak, adsorption of 
GNRs on Pd or Ti leads to the strong orbital mixing 
of the electronic states at the place of contact. The 
DFT calculations performed in this work show that 
the so-called metal-induced gap states appear in the 
energy gap of GNR at the interface between nanribbon 
and metal. These states can effectively penetrate on 
the large distance inside GNRs that can lead to the 
shortening of the metallic contacts and this effect limits 
the application of small GNRs in the devices. 

Presently GNRs can be prepared even with atomic 
precision in different ways [258-267 . Fig.|l7|]A) shows 
the STM images of (a) AGNRs and (b) chevron- 
type GNRs, respectively, synthesised on Au(lll) 
via bottom-up approach from different molecular 
STS measurements performed on 


the gap of 2.3—2.7 eV as deduced from the quasiparticle 
GW calculations corrected for the image charge in 
metals |248| . 

Later the combined ARPES and inverse photoe- 
mission (IPES) measurements were performed on these 
GNRs [ARPES data for 7-AGNRs: Fig. [T^B); IPES 
data: not shown] [248 269 . In this case the alignment 
of GNRs on the macroscopic scale is required and the 
stepped Au(788) surfaces were used for this purpose. 
This surface consists of the 3.83nm-wide {111} ter¬ 
races which can be used for the template growth of 
both types GNRs, straight ANRs and chevron-type. 
These experiments yield 2.8±0.4eV, 1.6±0.4eV, and 
3.1 ± 0.4 eV for 7-AGNRs, 13-AGNRs, and chevron- 
type GNRs, respectively. Modelling of these GNRs 
on Au(lll) within DFT and many-body electron ap¬ 
proaches gives a band gap of 2.85 eV and 2.96 eV for 


and electronic properties of nanoribbons and their 
dependence on the edge termination were studied 
via combination of STM and AFM with the CO- 
terminated scanning tips. It was found that initially 
synthesised GNRs have H-terminated edges and in this 
case the bonding to the Au(lll) substrate is weak 
that allows to move them with the STM tip. At 
the same time the STS measurements at the edge 
of GNR show the existence of the so-called vibronic 
tunnelling mode. In further experiments, the structure 
of GNR was modified via removing one of the H- 
atoms at one of the its end as was clearly identified 
via combined STM/AFM measurements - the bond 
formation between edge C-atom and the underlying 
Au was found. This modification leads to the blocking 
of the mobility of GNR on the surface and also to the 
drastic changes in the STS spectra. In this case the 
STS spectrum measured at the modified edge shows 
no fine structure (only very broad peak) and the one 
measured at the non-modified edge demonstrate the 
strong suppression of the vibronic mode. These results 
demonstrate the importance of the edge termination 
and the possible formation of the GNR-metal contacts 
for the modelling and interpretation of the transport 
properties of nanoribbons. 

There are two common ways to synthesise 
graphene nanoflakes (GNFs) or quantum dots (GQDs) 
on metallic surfaces: (i) the so-called temperature 
programmed growth (TPG) when on the first step the 
hydrocarbons (C 2 H 4 [60| 271 272 


or coronene 


60 


which are in rather good agreement with experimental 
values. These calculations indicate the charge transfer 


are predeposited on the metallic surface and then 
this system is annealed at elevated temperature or 
(ii) unfolding the predeposited Ceo molecules as 
demonstrated in Ref. [97 . The upper panel of Fig. 18 
shows 3D images of graphene flakes and quantum dots 
obtained in STM (a,c) and in combined AFM/KPFM 
(b) measurements. These GQDs were formed on 
Ir(lll) via TPG method and C 2 H 4 was used as 
a precursor. Combined AFM/KPFM measurements 
[Fig. 18 upper panel, (b)] yield the LCPD value 
of « 600 meV, which can be compared with the 
measured work function difference between graphene 
and Ir(lll) of 1.6eV and 1.1 ± 0.3eV obtained from 
LEEM 273] and STS 274 experiments, respectively. 
The observed discrepancy can be assigned to the 
smearing effect of the macroscopic scanning tip during 
KPFM measurements. This synthesis method of 
graphene nano-dots (GNDs) leads to the formation 
of well formed graphene islands with straight edges 
oriented along main crystallographic directions of 
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Figure 18. Upper panel: 3D representation of data for GNFs and GNDs on Ir(lll) obtained in (a) STM (180 x 180nm 2 ) and 
(b) combined AFM/KPFM (150 X 150nm 2 ) measurements. In (b) the topography of the system measured by AFM is overlaid by 
the respective KPFM signal measured simultaneously, (c) Atomically resolved STM image (20.2 X 20.2nm 2 ) of the single GQD 
on Ir(lll). Middle panel: (a,b,d) experimental and (c) theoretical STM images of the H-edge-terminated graphene nanoflakes on 
Au(lll). Bottom panel: (a) STM topography of an approximately 400 x 160nm 2 large flake used for dl/dV mappings, (b) dl/dV 
map of a large graphene flake marked by “b” in (a), (c) FFT images at selected bias voltages, (d) dl/dV map on the graphene flake 
with atomic resolution, (e) FFT maps of the image (d) at different bias voltages, (f) Plots of the Au(lll)- and graphene-related 
dispersions of the electronic states. The k values are plotted with respect to the T-point in the case of the surface state and with 
respect to the A-point in case of graphene. Data are reproduced from Ref. [202] with permission. 
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the Ir(lll) surface. The high quality of such 
islands allows the simultaneous atomically-resolved 
STM/AFM imaging of graphene and metallic surface 
[Fig. [l 8 | upper panel, (c)] giving a possibility to 
carefully trace the crystallographic and electronic 
structure of nanostructures on the atomic scale. 

The electronic structure of such GNDs on Ir(lll) 
was extensively studied by STS in several recent 
works 271,272 275 276 . These GQDs have exclusive 
zigzag edges. The presence of edge-states was not 
detected 277 as supported by DFT calculations and 
it was explained by a hybridisation of the GQDs p z 
orbitals and the substrate valence band states (here: 
Ir 5d z 2 surface state). Such interaction gradually 
decreases in strength from the edge towards the 
centre of the GQD. It is interesting to note, that 
although these works present the similar experimental 
observations, they do not provide a clear explanation 
of these effects (effects of quantization as well as 
the explanation for the extracted dispersion of the 
electronic states, E(k)) and did not give an answer 
about contributions of the electronic states of graphene 
and metallic substrate in the tunnelling current and, 
hence, in the observed effects. 

The properties of such GNDs can be tailored 
in different ways (intercalation, edge termination, 


manipulations, etc.). For example in Refs. 202 209 


the GQDs/Ir(lll) was intercalated either with Au or 
oxygen, respectively. In the first case, Ref. (202 


the thick layer of gold (50 — 100 A) was deposed on 
GQDs/Ir(lll) prepared by TPG and then system 
was annealed that leads to the formation of graphene 
flakes on Au(lll). The quality of this flakes is 
very high as can be depicted from Fig. 18 (middle 
and bottom panels). These STM images (middle 
panel, a) give a possibility to simultaneously resolve 
a herringbone structure of Au(lll), moire structure of 
the graphene/Au(lll) interface (as due to the lattice 
mismatch between two materials), atomic contrast of 
the graphene layer, and the edge-scattering effects in 
the GND. Comparison of the experimental STM data 
and the results obtained in DFT simulations (middle 
panel, b-c) demonstrates very good agreement for the 
H-terminated GNFs. The hydrogen termination and 
the weak interaction between graphene and Au(lll) 
are supported by the possibility to move flakes with the 
scanning tip as it was demonstrated in the experiment. 
Such possibility is absent for the graphene nanoflakes 
on Ir(lll) due to the C-Ir bonding at the edges as was 


demonstrated in Ref. [168 . 

The STS measurements performed on the 
GNFs/Au(lll)/Ir(lll) system allow to separate con¬ 
tributions in the tunneling current from graphene and 
from the substrate. Such experiments were performed 
on GNF, which is marked in the STM image shown in 


Fig. 18 (bottom panel, a) by letter “b”, and its STS 
map measured at the bias voltage of —75 mV is shown 
in (b). These data show the characteristic standing 
wave patterns and the corresponding ring structure 
in the FFT images obtained from the measured data 
at the different bias voltages (bottom panel, c) is un¬ 
ambiguously assigned to the surface state of Au(lll). 
The obtained dispersion of the electronic states, E(k), 
where qAu,gr/Au = 2fc, has a parabolic dispersion (bot¬ 
tom panel, f) with the effective mass to* = 0.26TO e . 

Analogous STS measurements performed with 
atomic resolution [Fig. [l8] (bottom panel, d-e)] show 
in the FFT images additional structure: (i) six spots 
corresponding to the reciprocal lattice of graphene, 
which are superposed by the spots originating from 
the reconstruction of Au(lll) and the moire structure 
(large hexagon); (ii) ring-like features which build a 
hexagon corresponding to the (V3x v / 3)R30° structure 
in the real space and related to the intervalley 
scattering (gi nter = FA' — 2k). These data obtained 
for different bias voltages allow to plot the dispersion 
of these electronic states [red squares in (f)]. The 
linear fit of these data gives the Fermi velocity of 
vp = (1.1 ± 0.2) ■ 10 6 m/s and the position of the 
Dirac point at E — Ep = 0.24 ± 0.05 eV, i. e. graphene 
is p-doped as found in the earlier photoemission 
studies of the graphene/Au(lll) system 159 . These 
measurements performed on the same graphene flake 
show the spectroscopic features from the surface state 
of Au(lll) and from the intervalley scattering of the 
graphene Dirac fermions allow to give an answer about 
different contributions in the imaging of graphene on 
metals [20^[20^[27l|[272|[275|[276l[278l . 

Later experiments performed on the GQDs/oxygen/ 
Ir(lll) system 209 give the Fermi velocity of vp = 
(0.96 ± 0.07) • lO^m/s and the position of the Dirac 
point at E — Ep = 0.64 ± 0.07 eV for graphene-related 
states and these values agree well with the data ob¬ 
tained in ARPES experiments in the same and similar 


works 163 209 


7. Conclusions 

The present manuscript reviews the recent progress 
in the studies of the structure and the electronic 
properties of the model graphene - metal systems. It 
is discussed that the properties of graphene are defined 
by the interplay of morphology of the system and the 
overlap of the electronic valence states of graphene 
and metal at the interface. These effects define the 
modifications of the dispersion of the electronic states 
of graphene around the Fermi level, like the shift of 
the Dirac cone (doping), destruction of the Dirac cone 
due to hybridization of the valence band states of 
metal and graphene, or/and band gap opening in the 
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electronic spectrum of the graphene ir states due to 
the violation of the sublattice symmetry in graphene. 
Considering all these effects one can conclude that 
electronic spectrum of graphene, which is in contact 
with metal, is always strongly disturbed leading 
to the loss by graphene its unique properties that 
influence the transport properties of the devices where 
graphene/metal junctions are present. Here we try to 
summarise some ideas that might help to overcome 
these difficulties via application of different methods 
that tailor the properties of the graphene/metal 
interfaces or via formation of nano-objects on the basis 
of graphene on metals. 
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